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THE  THEORY  OF  SURFACE  WAVE  DIFFRACTION  BY 


SYMMETRIC  CRUSTAL  DISCONTINUITIES 
By 

J.  Kane 
and 

J.  Spence 
summary 

A  major  barrier  in  comparing  seismic  theory  with  observed  wave 
trains,  stems  from  the  fact  that  elastic  wave  characteristics  are  in¬ 
fluenced  by  discontinuities  along  the  propagation  path,  and  any 
understanding  of  such  signal  corruption  would  require  a  knowledge  of 
the  diffracted  fields  at  the  appropriate  obstacles.  However,  even  a 
relatively  simple  crustal  feature  such  as  a  discontinuous  change  in 
terrain  presents  major  difficulties  if  the  relevant  problem  of 
diffraction  in  a  wedge-shaped  region  is  considered.  Although  some 
first-order  calculations  have  been  made  by  Lapwood  (1961),  Kane  and 
Spence  (1963),  and  Hudson  and  Knopoff  (1963).  the  theoretical  discussion 
of  the  diffraction  effects  are  hampered  by  the  intractability  of  the 
associated  boundary  value  problems.  In  this  report,  we  show  how  one 
can  take  advantage  of  symmetry  considerations  and  variational  techni¬ 
ques  to  rapidly  estimate  reflection,  transmission,  and  conversion 
coefficients  for  elastic  wave  diffraction  at  syrisot r ;e  wedge-shaped 
obstacles.  In  Hart  I,  we  illustrate  the  Ideas  by  a  discussion  of  the 
vector  problem  of  Rayleigh  wave  propagation  along  the  faces  of  an  clastic 
wedge  with  tree  boundaries.  In  Part  II,  we  analyze  the  scalar  problem  of 
multi-mode  Love  wave  diffraction  in  i  symmetric  layered  wedge. 
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FART  1.  RAYLEluli  WAVES  ON  AN  ELASTIC  WEIXjE 
i  .  Fun  camenta  1  equal  inns 

The  tremors  slu.v.w)  of  an  elastic  sol  ic  cnaracter lzed  by  the 
Lame  parameters  •  ,  n,  and  density  •"  .  can  be  dei  ived  from  a  scalar 
potential  rU.v.z),  and  a  vector  potential  ,  *  U.y.z)  , 

(x , y , z ) ]  by  the  relation 

s(u,v,w)3V,c«-V-'lr.  II. I) 


For  two-dimensional  motions  which  are  independent  of  the  /.-coordinate, 

both  t  and  «r  are  but  functions  of  x  and  y,  or  ?  *  :U.y),  and  V  *  . 

Furthermore,  we  can  neglect  pure  distortions  by  setting  .  *  •  J  0, 

A  > 

so  that  the  vector  potential  *•'  *  f  [U.  ,s 

characterized  by  one  scalar  cwnponent  and  the  Subscript  on  i  can  be 
dropped  without  contusion.  If  we  assume  that  the  vibrations  are 
harmonic,  we  can  suppress  a  t  i~*e  factor  e  .  .«nd  r.  ran  be  shown  th.it 
*,  and  the  ? -component  of  the  vector  potential,  sal  is: v  the  reduced 
wave  equations 


IV 


2 


O  7  *  O. 
c 


k" 

c 


t  *u ) 


,1.21 


(V 
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k‘)  t  * 
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Once  the  potentials  :  and  *  arc  known,  the  displacement  vector  s  is 
given  by  (1.1).  .and  the  resultant  stress  dyadic  (7,?>  can  be  tivun 
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in  symbolic  notation  as 

-  C  ,0  "  '  ^  V-s  *  e  (.Vs  ♦  sV) 


u  .4) 


where  3  is  Lhe  unity  dyadic,  the  ideinfactor. 

d .  Die  Havleiih  wave  potentials 

A  time  harmonic  Rayleigh  wave,  or  fi-wave  i  ai  brevity,  is  com¬ 
prised  of  a  pair  oi  exponential  solutions  o:  (1.2)  and  (1.3), 

«  i(*  ,  -j  ) .  If  the  elastic  solid  lies  within  the  halt-space 
K  it 

y  s  0  (cf .  Figure  la),  then  these  solutions,  r.  polar  coordinates 
x  ■»  r  cos  u,  y  *  r  sir.  0,  assume  tlio  torn 
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•  (tij- 


tl.  a) 


t  .  .  0  ) 


whe  re  li  t  ill.  o  ’■  p 


v  t!  lit.  :>  :s  ,;!Vc! 


"Hill  ,1  -(V  v  >  ‘  :i  o  ■  1  e  . 


"01  i -  1  v  I  ‘  i 


,,m!  "  ,  t  he  ■  !.-n  ;  !  oCi  !  '  hi 


I  the  >h«.  ar  , 
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and  is  a  dimensionless  distance  parameter 


v  ■  — 

**  K 


(I  .  HI) 


The  parameters  v  and  v  represent  the  velocities  of  the  compres- 
sional  and  shear  body  waves  respectively,  foi  a  oven  Poisson's  ratio 


2(  •  ‘  a) 


f  1.11) 


one  needs  choose  the  hayloigh  wave  velocity  vR,  wlncti  is  less  than  v^, 
so  that  the  stresses  induced  by  and  tR  vanish  at  the  surface  y  t' 


(  •  .  *  )  •*  u , 

H  V 


y  -  0. 


(1.12) 


With  this  choice  o:  vR.  the  CR  and  iK  it i veil  bv  ll.1))  and  tl.6)  are  the 
vector  and  scalar  potentials  character  i  Mint  a  Kayloinh  wave  travel  in* 
to  the  ih;Iu  with  unit  amp  l  Unde  .  (wC  shall  s|H;  ik  •’(  the  c  h:I  f  l  c  lenl 
of  the  compress  i  ona  1  potential  as  the  amplitude!.  ! '.  is  very  con¬ 
venient  to  note  that  we  call  reverse  the  direction  .  any  li.irnoni  Wave 
bv  the  operation  of  complex  cwii’.nMtlcir,  thus 


(  * 

I  L  «  .  U 


(1.1  II 


••  "S.'V*  <t 


*  -  1 0  ) 
H  <K  -  *»K  e 


t  I  .  I.’.) 


represents  a  Rayleigh  wave  travel  me,  in  the  <>ppo.-  a  reel  i  >:  v  th 

ampl  it  uric  R. 


3.  Formulation  of  the  boundary  value  problems 


a.  The  major  problem 

The  conundrum  posed  by  Kuvleigli  wave  dll  fraction  in  a  wedge- 
shaped  region  is  to  find  additional  solid  ions  / ^  and  v^,  of  (1.2)  and 
(1.3)  which  represent  d  1 1 1  ratted  fields  m  miciiur  of  the  wedge 
such  that  the  stress  dyadic  -(“^  *  r^.  >r  *  :(!'  vanishes  on  both 
faces  of  the  wedge  for  8-wave  excitation  along  one  face.  In  our  geun- 
etry,  (£f.  Figure  l.i)  ,  the  8 -wave  is  incident  from  infinity  along  the 
negative  x-axis.  We  shall  he  principally  concerned  with  calculating 
the  complex  amplitudes  of  the  reflected  and  transmitted  »-uaves  as  a 
function  of  the  wedge  angle  0  and  Poisson's  rate  This  task  .. 
be  substantially  cased  by  reducing  the  major  ;>t  hie-  to  a  pair  of 
minor  problems  involving  even  and  odd  syiirr.ctr  its .  (  ns  ><  i  the  in¬ 

cident  d-wave  along  the  left  wedge  face  to  be  me  sum  of  two  waves, 
each  of  half  amplitude,  likewise,  the  absence  of  anv  excitation  along 
the  right  wedge  face  is  equivalent  to  a  pair  oi  u  ident  '"--waves  along 
ll,  each  of  half  amplitude.  Put  opnosed  in  Sj  .  .  icf.  Figure  Mil.  Hie 
four  partial  R-waves  can  be  separated  into  two  groups:  First,  a  pair 
of  r-vavos  on  either  taco  of  like  parity  whi  'h  serves  a  u1  excital i  ai  f 
what  wa  call  the  even  uimr  problem,.  Second .  another  p..  i  r  >i  1  wives  whose 
amplitudes  are  of  unlike  parity  which  constitutes  ih«  excitation  oi  the 
odd  minor  problem,  it  we  designate  the  diffru-’ ed  '  nhh  >■:  the  even  md 
odd  problems  by  the  subscripts  c  and  o  respectively,  then  since  we  are 
dealing  with  linear  cm- ml  ions,  the  desired  major  potentials 


and 
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*  can  be  expressed  as  a  superposition  of  the  minor  potentials 
d 


<1  we  "o  ' 


(1.15) 


A  =  iV 

*d  ‘e 


(1.16) 


and  likewise  the  displacement  vector  a  =  can  be  decomposed 

into  even  and  odd  components. 


b.  The  even  minor  problem 

In  the  even  problem,  the  wedge  will  suffer  only  even  dis- 

★  — 

placements  s  about  the  plane  of  symmetry,  and  as  a  result,  there  can 

-  e 

be  no  component  of  normal  displacement  along  the  plane  of  symmetry 

atO-6V  2  -  n.  It  follows  that  the  even  problem  is  equivalent  to 

finding  the  potentials  in  a  bisected  wedge  with  one  face  free  of 

stresses  which  supports  the  incident  Rayleigh  wave,  and  the  other  lace 

so  constrained  that  the  norma  1  displacement  vanishes  there.  Thai  is, 

we  seek  solutions  .*  and  v  of  (1.2)  and  (1.3)  in  a  wedge  of  hnif- 
c  e 

angle  6^/2 


EVEN 


trK  ’  *e*  ‘r  ’  re>  ‘ 


0  »  ii 


-  —  (:  *  :  )♦  —  (t,  *  v  )  *  2  -  a 

r  ao  R  o'  or  R  c 


(1.17) 


(l.lh) 


c.  The  odd  minor  problem 

by  the  same  arguments,  the  odd  problem  which  involves  :  , 

is  a  complement arv  version  in  the  bisected  wedge,  wherein  the  tangent iai 

*N.B.  The  displacements,  and  the  compress ional  potential  will  he  even 
about  the  plane  of  symmetry,  but  the  shear  potential  will  he  an  odd 

function,  and  v  i  .'c  versa  for  the  odd  problem. 
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d isplaccmencs  most  vanish  identically  along  the  plane  of  symmetry,  i.e.. 
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d .  rh v  He i  1  c  cl  i  vi.  Coe  1 1  i  c  lens  s 

For  either  the  even  or  'lie  add  j> r  >b  1  e  .  'he  .»olut  mu  will  contain 
a  roilected  Itaylcw.h  wave.  Let  _ l  &  1,  '!  and  (('9  7,  ')  he'  ihe  c on  - 

plc.\  ratio-  ot  the  retie  •  ec:  tu  the  incident  Rayleigh  wave  amplitude 
for  the  even  ami  odd  -nnoi  problems  in  the  hi. -verted  wedge .  The  re¬ 
flection  coefficient  Kt/y,  )  :or  the  original  'aj  r  p rub  lens  will  be 

i-t  t  • "  L  ■’ . ")  ’  ((v  ..oi.  .  (1.71) 


and  1  Kc  w  l  ra  i  he  e>v<-al.  t  r  an  smuts  tun  -k  1 1  i  icni  l  i^J) .  will  he 


lt^.-i  "  ;  t  (<'• 
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,1  i 


( y  . .  j . 


(  !  .  .»») 


For™.i  i  I  .1  Wit*  a:,  be-  V -1  it  •  ••  ••  e  at  F  i  gtirc  lb 

which  indicates  •  hat  ’he  -  oral  ret  !•  -  ‘  ;  1  1  4 ’’  **  rese.1- 

f  roc.  a  s-atu- 1  ;  Sit  ion  ot  the  partial  retie  ti  ti  *  el  1 1 .  lent  s  t(.:.(/  S  )  . 


and  the  l  raiisii. .  i  on  rot  if  icteic  .  r  r.. 


laeif  i nt e r l  v :  once  .1.'  *  .  (  ) 


/, .  The  Variati  oi.  il  i  •  i  t»le 
a .  I)  l  sens  si  i  w'l 

V.iriat  i.in  u  pr  -  i  ,r<  -  consist  of  a  •  sn-ing  i  suitable  trial 
fund  lot.  :o;it  a  in  in  .  u:  pec  I '.  led  e .  e  f  t  i  c  lev '  s  .  u.i  t  hen  chons  In;;  these 


-8- 


parameters  to  minimize  certain  quantities.  One  major  advantage  of 
the  variational  method  is  that  first-order  accuracy  in  the  trial 
function  usually  gives  results  which  are  accurate  to  second-order, 
because  of  the  stationary  character  of  the  approximation. 


A  natural  aperture  in  the  present  problem  is  the  plane  of 
symmetry  and  we  can  assume  it  to  be  illuminated  by  an  incident 
8-wave,  and  a  reflected  one  with  an  adjustable  amplitude.  In  the 
even  problem,  the  net  angular  displacement  must  vanish  along  the  plane 
of  symmetry.  A  unit  ft -wave  traveling  to  the  right  gives  rise  to  the 
angular  component  of  the  displacement 


inc 


,  r  8  '  i 

^  ^  oO  C  re  ■*’ 


(1.23) 


and  likewise  an  •’ -wave  of  amplitude^  moving  to  the  left  generates  the 
disturbance 


ref 

c  S0 


(s 

e  u 


inc . 


(1.2M 


which  apart  from  an  amplitude  factor  is  the  ur.picx  conjugate  of  (1.23). 
Only  If  there  is  no  discontinuity,  or  <9  it  can  we  make  the  angular 
displacement  of  the  trial  function 


inc 


,  inc,’ 
r  ( sn  ) 
e  0 


(1.25) 


vanish  for  all  r  along  the  plane  of  symmetry  bv  prnpcrlv  choosing 

4c  1* 

Otherwise  r,  i  a* ,  j  *  ,  and  no  choice  of  rc  can  make  sQ 
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vanish  at  more  than  an  isolated  set  of  points.  There  are  at  least  two 
ways  by  which  we  can  improve  matters:  We  could  use  a  more  complex 
trial  function  which  acknowledges  body-wave  contributions  to  the  dif¬ 
fracted  field,  or,  since  the  residual  displacement  sQ  is  explicitly 
known,  we  can  use  it  as  the  aperture  illumination  of  a  Green's  theorem 
type  calculation  to  correct  the  variational  es'imate. 

However,  the  practical  seismic  interest  is  in  the  realm  of  small  discon¬ 
tinuities,  and  for  this  case  we  shall  sec  that  the  elementary  trial 
function  yields  satisfactory  results. 

b.  Definition  of  the  scalar  product 

While  there  are  many  guages  by  which  s^.  can  be  minimized,  we  shall 

T 

choose  a  c  which  minimizes  s  in  the  mean  square  sense.  For  this  pur- 
c  u 

pose,  let  us  define  the  complex  scalar  product  of  two  functions  ulr.O) 
and  v(r,0)  to  be 

( u , v )  s  u(r ,0)  v  Cr.O)  dr,  ^  • 

O 

where  the  integration  is  to  be  carried  out  along  the  plane  of  symmetry 
0  .  if*.  1 2  —  n.  To  each  complex  function  ulr.O)  we  can  attach  a  positive 
definite  number,  the  norm  of  u  or  ||  u  |j  which  is  defined  to  be  Au  7u) . 
The  norm  ||  u  j|  dejxinds  on  the  wedge  angle,  and  is  to  b"  list  mguished 
from  u2  -  (u,u*)  which  Is  in  general  complex. 


c.  The  Even  Subsidiary  Reflection  Coefficient 

With  this  notation,  the  mean  square  value  of  the  angular  dis- 

T  . 

placement  of  the  trial  Lunction  sQ  is 


,  me, v  inc 
e  S0  )  ’  S0 


,  tnc  ’ 
( s  _  ) 

e  0 


26) 


Kn2=Gr 


a. 27) 


-  10- 


and  tVis  \  i  i  1  be  .1  .  uw.  1.  it  u.ij  only  it  Is  eheaeii  .. 


.  (&nn  «  — - 
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*  I  i 


( I .  2‘i) 


or  explicitly  m  terms  et  <(u),  9(0)  411U 
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Figure  2:  The  diffraction  coefficients  K  and  T  for  a  trial  function  con¬ 
sisting  of  but  an  incident  and  reflected  wave.  In  this  case, 
Poisson's  ratio 


1/4. 


ARG  R (CIRCLES)  ARG  T(CIRCLES) 
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With  an  explicit  .'g  and  rQ  at  our  disposal,  we  can  evaluate 
the  R(^,-)  and  T(^,")  germane  to  our  elementary  trial  function,  r.ir 
reference,  their  complex  variation  is  plotted  as  a  lunation  of  the  dis¬ 
continuity  angle*?  -n,  in  Figure  2  for  Poisson’s  ratio  ■  -  1/4.  The 
present  magnitudes  |r|  are  somewhat  smaller  .■>  ’-ha'  iven  by  earlier 
first-order  calculations  (.Kane  and  Spence  l1;'  >  .  which  do  not  simul¬ 
taneously  yield  K  and  T.  Since  we  evaluate  both  R  ntd  T  together,  we 
must,  in  effect,  withdraw  some  energy  from  the  iiaiisn.it ted  field  to 
allow  for  the  reflected  wave.  Furthermore.  is  the  nature  of  the 
variational  technique  to  underestimate  the  subsidiary  diffraction 

coefficients  ;  and  ;  since  it  only  yields  their  projection  in  the 
e  o 

sub-space  spanned  by  the  trial  function. 

Although  the  analysis  is  certainly  valid  for  a  small  enough 

discontinuity  in  wedge  angle,  the  utility  of  the  procedure  can  mu  be 

established  until  there  is  swne  estimate  o:  the  errors  committed.  A 

feature  of  the  present  procedure  is  that  it  suggests  a  natural  -.wage 

*1  T 

of  the  accuracy.  While  and  :q  are  so  chosen  that  ||  s()  ||  and  !|  sr  || 
arc  minimized,  both  s^  and  s^  are  non-zero  along  the  aperture  plane. 
These  residuals,  which  are  explicitly  known .  can  not  be  laitncr  re¬ 
duced  without  introducing  new  features  such  as  body  wu  -mr. but  ions 
into  the  analysis,  since  ll-|R| 2-|t| 2]  represents  in..'  fra  i  ion  of 
energy  unaccounted  for.  we  can  estimate  the  need  for  i.npw.ng  the 
calculations  by  examining  this  quantity* 
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This  error  estimate  is  a  very  generous  one  because  only  part  of  it 
implies  higher-order  corrections  to  R  and  T,  the  remainder  representing 
energy  which  is  accounted  for  by  R-wave  to  body  wave  conversion.  The 
data  of  Figure  2  shows  that  if  |^_n  |  ‘  10  ,  the  present  analysis 
accounts  for  at  least  92  percent  of  the  energy,  and  therefore  the  theory 
should  not  require  further  improvement  within  this  range.  We  can  also 
compare  the  present  theory  with  experiment,  but  we  must  be  very  careful 
if  we  do  so,  because  there  are  fundamental  distinctions  between  analysis 
in  the  harmonic  domain  and  pulse  measurements  i£f.  Appendix). 
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PART  II:  LOVE  WAVES  ON  AN  ELASTIC  WEUGE 
1.  Introduction 

A  layered  solid  can  support  surface  waves  which  .uc-  •'  <  >  >n"  i 
polarized  shear  waves  trapped  in  the  superfici  il  layer* 

Since  these  Love  waves,  as  they  are  known,  have  i  impress  tonal  com¬ 
ponent,  it  is  not  necessary  to  introduce  potent  ial.!,  and  it  is  possible 
to  work  directly  with  one  scalar  function  w(x,y),  the  /-component  of 
the  displacement  vector 


s  •  [0,  0,  w(x,y) ] #  (2.1) 

Within  the  E^  layer,  w  satisfies  the  wave  equation 

(V2  ♦  kp  w(x,y)  ■  U,  kj  *  1  ■  j/Wj.  (2.2) 


and  within  the  E^  substrate,  w  obeys 
(V“  ♦  kp  w(x.y)  *  t>. 


k2  14  ‘  "2/u2‘  (2.1) 


At  the  free  surface  of  Ej ,  the  stress  dyadic  6( s') must  vanish,  which 
will  he  true  provided  that  the  normal  derivative 


dw 

On 


0 


(2.4) 


vanishes  there. 

We  assume  the  E^E^  Interface  to  be  welded  so  that  the  displacement 
and  normal  stress  must  be  continuous  across  this  boundary.  Along  wedge 
face  A,  these  conditions  will  be  satisfied  provided  that  the  Love 
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waves ,  or 


•fj-waves  of  amplitude  A.  have  the  form 


n — t  'V 

A .  cosvk,  -  vye  , 


'I'  i* 


0  y  -  H 


A.fl  ■ 

t  i 


2-  ■* 


1  • .X 
1 


A.  cosVk .  -  ‘..II  exjftX  -  k_y  e  ,  y  -  -  11, 

l  111* 


(2.6) 


(2.6) 


and  the  propagation  constants  •  are  the  real  r  ts  of  the  ;>eriod  eiiuation 


(2.7) 


If  the  shear  wave  is  to  be  trapped  in  the  layer,  or  if  the  Love  wave  is  to 
propagate,  we  need  ] k 2 1  •-  |'.|  '  |k(|.  For  any  thickness  however  small, 
there  is  at  least  one  root  j  corresponding  to  an  acceptable  solution  — 
i he  fundamental  Pj-wave.  As  the  acoustic  thickness  k^!  increases,  other 
modes  can  propagate.  In  our  discussion,  we  shjll  assume  that  the  thickness 
is  such  that  only  two  modes  propagate,  the  fundamental,  and  oik*  harmonic: 
the  P. -.  wave  with  a  propagation  constant  ' Hie  analysis  proceeds  in  a 
similar  fashion  if  an  arbitrary  number  of  modes  can  propagate. 

2.  Formulation  of  the  boundary  value  problem 

We  assvime  that  an  Pj-  wave  is  Incident  along  oik  Jute  oi  a  syirmet  r  leal  1  y 
layered  wedge.  At  the  discontinuity,  four  surface  v.jv.  will  be  excited: 

A  reflected  and  transmitted  i*  -  wave  with  amplitude  cuciiicients  and 
respectively,  and  reflected  and  transmitted  waves  whose  amplitudes  are 
the  conversion  coefficients  R  and  respect ive ly .  Our  task  will  be 
to  determine  these  diffraction  coefficients  as  functions  of  the  wedge 
angle ^  the  layer  thickness  h.  and  the  elastic  constants  s  ,  u  ,  n2> 
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By  the  same  argument  as  in  Fart  I,  we  can  add  and  subtract  a  symmetric 
r^-excitat ion  on  the  right  wedge  face  which  leads  us  to  consider  a  pair  of 
even  and  odd  problems  in  a  bisected  wedge,  since  Love  wave  diffraction 
is  a  scalar  problem,  the  subsidiary  boundary  conditions  along  the  aperture 
or  plane  of  symmetry  are  simply 


EVEN : 


ODD: 


1  dw 

- -  -  u . 

r  dO 


u  -  (9/2  -  n  (2.8) 

0  »  0/2  -  n  (2.9) 


for  the  even  and  odd  problems.  The  trial  function  will  consist  of  an 

*  * 

incident  J^-wave,  and  reflected  and  1,,-waves,  with  unknown  ampli- 

tudes.  If  we  denote  the  subsidiary  reflection  and  conversion  coefficient 

for  the  even  problem  In  the  bisected  wedge  as  r^  and  r  p ,  and  similarly 
0  0 

r^  and  for  the  odd  problem,  then  the  desired  maj  >r  cceff iclents 
are 


Rn  '  >  uu  *  rn\ 


(2.10) 


R12  *  *  U12 


0  , 

rl2  ) 


(2.11) 


'll 


i  ,  e  0 
V  (rn  -  r11)> 


(2.12) 


T  I/O  0  . 

T12  *  **  r12  *  r  1 2  * 


(2.13) 


As  in  Part  I,  we  shall  detenu ine  these  coefficients  by  a  variational  pro¬ 
cedure  which  ignores  body  wave  contributions. 
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3.  Sol ut ion 


In  the  odd  problem  we  shal 1  choose  and  r^y  so  that  the  residual  variation 


•:(r)  alont;  the  aperture  plane 


e(r)  -  0  ♦  r°  »"  *  r“  «! , 

etfi  1  rH  I  *12  1’ 


V  *  &  /  2  -  It  (2.14) 


is  as  small  as  possible  in  the  mean  square 


;  the  same 


definition  of  scalar  product  as  in  l’art  l,  we  h  ivi 


Or)  II2  -  (.,  .  r“,  *,  •  »;  •  rf,  «*2  ),  12.15) 


o  o  ,  r 

and  this  expression  will  be  a  minimum  if  and  only  li  rH  and  r12  saiisty 
the  normal  equations 


(YV  *  rn  lY  V  *  r!2  °Y  V  *  U’ 


(2.16) 


(p  i»  )  ♦  r°  (V'*  )  •  r°  (i*  .  1*  )  «  U. 

'  1 ’  2'  11  l’  -  12  22 


(2.17) 


Squat  ions  (2.16)  and  (2.17)  can  immediately  be  •  ved  :'oi  r^  and  r^ 


■  >V  V  ‘V  V 


(2. IS) 


<Y  V  lV  V  • 


(P  .  I*  )  (\>  V  ) 
i  '  i  •  r  i  •  i 


l  2. 19) 


(i-i .  v,)  oy  V  • 
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DET 


(Q  a  )  to  >'  ) 

v  1  ’  12’  ! 

it  if 

to  v*  )  (V  .  i‘  ) 

v  j  .  2  v  2  2 


(2. 20) 


e  o 

In  the  same  lash  ion,  lor  the  even  problem  we  .1  el. nose  r  ^  ^  and  r,9  as 


1 1  DET° 


12 


!!i 

o'.1  _ 

1  __T 

l 

-*2 

dO  ’ 

r  dO  . 
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and 


c  .  _  I 

V  c 

1  DETC 
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,  ovl, 

1 _ 1. 
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J.  _ 1 
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whe  re 


t)ETC 
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2i 
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OG  ’ 
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.10  / 
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00 

Huh  this  knowledge,  the  reflection,  transmission  imi  r» son  eoc£ £ ictenis 


are  given  by  (2. lit)  -  (2.13). 


C.  Discussion  of  the  results 


1 .  Numerical  Data 

We  have  used  the  pro  tx'din  •  i  omul  as  i  'a  leu  l  ue  the  diffraction 
coefficients  lor  .in  E^-laver  am!  b  , -Substrate  for  which  k  /k.j=  1.297 
and  ^  2.lc9.  The  phase-  *nd  gr  up  vel  >ei*  ten  for  this  case  have 

been  given  by  stonely  and  are  available  in  i  •>  aid  reference  Ip. 213, 

Ewing,  Jardetskv.  and  1'ress).  figures  4  through  illustrate  the  varia¬ 
tion  of  the  magnitude  f  tin  diffraction  cuef  f  ieient  s  which  are  even 
functions  of  the  discontinuity  angle  -n.  Hie  curves  are  indexed  by 

four  values  of  the  d mens i on  less  parameter  jit,  namely  i,6,7,b;  if 


•l" 

-  3,  then  the 

sc  cone1 

mode  : 

i  s  just  above 

cut  -  t  !  ,  and  if  •  jll 

-  8. 

the 

third,  mode  is 

Just  be 

l  ‘  V  e*\ 

,t  -  :  i  .  It  i. 

.  veiv  interesting  to 

note 

that 

the  c Oliver.',  . 

on  c  ». :  1 

i  :  tent 

! | j  v\.eed- 

the  rot  le  t  ion  oet  f 

ic  tent 

that 

is,  there  i» 

l  l elide* 

r.  v  ! 

•  •  he  ellergv 

ti  "itt  inue  •  propa 

■ate  in 

the  Sane  direction  <-v:  I  *  uc  e’aSitate  .  t  •  a’  ter  :  >'.i!  charac- 

t  e  r  i  s  1 1  c  s . 

2.  Interpret  at i on 

If  we  eauparc  any  tv  wive  <  :  idem  i--.il  a  •  '  e  i  ist  ies.  then 
their  relative  energy  is  prop,  rt  i  nil  •  :»e  e  %  •!•  are  of  any 

correspond  in,  .i-r.i ,--l  i  • .  ude  .  1  •'  •  ‘v  i  hand  .  a.;  a  t  -pare  the 

energy  in  the  1  uiida-sent  a :  X  wav  •  ■  hi!  of  an  2,-vu  *  iirst  har¬ 
monic,  we  need  make  s  « *  ’  irtne;  -alcul.it  ions.  Wi  *  •  ■  In,,  m  general  it 

let  us  s  IK:  i  a  1  t  -  li  •  s  >•  li  the  her  i  ?  ail  i  -  f  ace  lor  which 

the  X  -waves  are  giver  <  1  Hlv  bv  (2.5)  and  <  2.6)  .  and  evaluate  the 

i 

scalar  product  a  inn.  ''a  wavefront  v  *  li.  Thus  |  \  |  '  ||  X  ||  “  represents 

1  » 
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Figure  4:  The  magnitude  of  the  Love  wave  diffractiur  coefficients  for 

>.jH  -  5.;  they  arc  even  functions  of  the  discontinuity  angle 

•  \ 

O-  "•  In  this  case  k^kj  "  1.297,  M^/Mj  “  2.159;  the  normal¬ 
ization  value  ,\  -  2.U. 

Figure  5:  The  magnitude  of  the  Love  wave  diffraction  coefficients  for 
M  "  6.;  they  arc  even  functions  of  the  discontinuity  angle 

n.  In  this  case  k,,/kj  -  1.297,  "  2.159;  the  normal- 

izat Ion  value  N  -  1.2. 

Figure  6:  The  magnitude  of  the  Love  wave  diffraction  coefficients  for 
• jll  '  7.;  they  are  eve-  functions  of  the  discontinuity  angle 
y *  n.  In  this  case  k^  k^  *  1.297,  u,,/a  -  2.159;  the  normal- 

i zat ion  value  S  »  l . 1 . 

Figure  7:  The  magnitude  of  the  Love  wave  diffraction  cooff icients  for 
*jll  »  8.;  thev  are  even  functions  of  the  discontinuity  angle 
(£)  -  n.  In  this  case  k^/ k j  *  1.297,  i^/.j  *  2.159;  tit  normal- 
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the  mean  square  energy  flux  transported  by  an  J  .-wave  of  amplitude  A,. 

If  we  denote  the  group  velocity  of  an  I, -wave  as  t  11  follows  that  the 

rat  lo 

'l  1  fl  X  (2.24) 

I  A,  |  2  II  t£  ||2 

compares  the  power  flow  ot  an  A.Tj-wavc  to  an  r>  >j  f - wave .  In  particular, 

a  mode  near  cut-off  behaves  like  an  unbounded  plane  wave  in  the  b ^ - me d i um , 
hence  such  a  wjve  can  carry  large  amounts  of  power  even  if  its  amplitude 
is  deceptively  small.  As  a  result,  if  vc  arc  to  discuss  power  transfer, 
we  should  renormalize  the  amplitudes  of  the  conversion  coefficients 


NK 


12* 


V 

r'.2 


NT 


12 


(2.25) 


where 


(2.26) 


so  that  |KJj2l 2  and  1^,1 2  ire  proportional  t<  '  >e  power  transferred  bv 
the  diffraction  of  a  Tj-wavc  of  unit  amplitude  i  wedge  discontinuity. 
The  appropriate  values  of  N  tor  the  previous  numerical  example  arc 
cited  In  the  captions  of  figures  4-7. 


In  a  fashion  similar  to  the  tint  analysis  ot  I’.irt  l.  ihe  function 
[l.|K  |2.|R^2j2-|T1,|2-|T;i*2|2]  represents  the  amount  I  ambiguous 

energy.  These  values  are  more  satisfactory  in  the  present  analysis 
than  in  fart  I;  this  can  be  explained  by  the  fact  that  we  have  a  more 
flexible  trial  fun. lion  nice  we  can  vary  the  coelf icicius  of  two 


reflected  modes . 


APPENDIX 


Love  wave  diffraction  coefficients  would  be  very  difficult  to 
measure  in  the  laboratory,  but  the  techniques  of  two-dimensional  model 
seismology  offer  a  means  of  determining  Rayleigh  wave  reflection  and 

transmission  coefficients  with  an  accuracy  of  about  10-20  per  cent. 

^  v  Y 

The  present  theory  and  experiment  agree  if  -  «,  but  outside  this 
range,  there  are  experimental  features  which  are  not  duplicated  by  the 
results  of  the  present  elementary  variational  procedure.  The  analy¬ 
sis  in  the  harmonic  domain  could  be  refined  by  employing  various 
devices  to  reduce  the  amount  of  unexplained  energy.  Such  calculations 
would  probably  require  substantial  effort,  anti  the  idealized  formula¬ 
tion  of  the  present  problem  should  be  reviewed  if  the  labor  is  to  have 
relevance  to  puise  measurements. 

There  are  major  distinctions  between  analysis  in  the  harmonic 
and  time  domain.  For  example,  whereas  a  harmonic  Rayleigh  wave  is  a 
uniquely  defined  entity.  Friedlander  (1%8)  has  pointed  out  that  a 
Rayleigh  pulse  can  assume  a  variety  of  waveforms.  Furthermore,  any 
Rayleigh  pulse  can  not  have  a  sharply  defined  wavefront,  and  theoreti¬ 
cally  must  give  infinite  advance  notice  of  its  arrival,  unless  it 
merges  continuously  with  a  precursor,  typically  the  shear  psuedo-surface 
wave  (Cagniard  1939).  Although  the  amplitude  of  this  shear  wave  decays 
Vlth  distance,  its  integrated  flux  remains  constant.  It  is  difficult 
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Figure  8:  The  diffraction  coefficients  R  and  T  for  a  trial  function 

consisting  of  an  incident  and  reflected  Rayleigh  wave.  How¬ 
ever,  the  shear  coefficient  7  of  the  incident  Rayleigh  wave 
has  been  incremented  by  a  factor  (1  +  cos  +  cos  ’JS ) .  In 
this  case,  i’oisson's  ratio  n  =  1/4. 
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to  separate  the  far-field  effects  due  to  the  arrival  of  a  Rayleigh  pulse 
and  its  shear  companion  at  the  second  wedge  face.  In  other  words,  in 
addition  to  Rayleigh/Rayleigh  interactions,  there  will  be  some  shear /Rayle igh 
conversions.  What  contributions  might  this  shear  wave  introduce?  Whereas 
we  can  not  give  a  rigorous  answer  to  this  we  can  however  make  a 

rough,  but  simple,  estimate. 

We  first  note  that  if  the  wedge  angle  9  is  n  or  11/ 2,  then  we 
would  expect  little  or  no  shear/Rayleigh  conversion.  In  the  first  case, 
there  is  no  discontinuity,  and  the  second  case  corresponds  to  a  geometry 
for  which  the  shear  wave  is  essentially  normal  to  the  second  wedge  face, 
and  we  know  that  for  normal  incidence,  a  shear  wave  is  reflected  as  a 
shear  wave.  Then,  from  Equation  11.6),  we  note  that  we  can  easily  add 
some  additional  shear  potential  to  the  original  excitation  by  incrementing 
the  Rayleigh  wave's  shear  coefficient  by  an  additional  contribution 
f<©»  depending  on  the  wedge  angle 

*  •  [  1  •  f  (f*)3 

Quite  arbitrarily,  we  have  chosen 

-  <-*  IT-. 

(  {+*'.)  “  cos*"*  +  cos  * 

whose  sole  merit  is  that  it  is  the  simplest  function  we  could  think  of 
that  vanishes  for&y-  «,  "/2.  It  is  then  a  trivial  matter  to  repeat 
the  calculations  appropriate  for  Figure  2,  and  the  results  are  plotted 
in  Figure  8.  The  shaded  area  indicates  the  range  of  exper imentnl  points 
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as  measured  by  Knopoff  and  Cansji ( 19 6Ci  ,  deBremaecker ( 1958)  .  and 
Viktaruv,  (1958).  Of  course  there  is  limited  justification  for  this 
heuristic  procedure,  but  it  is  remarkable  that  with  this  naive  device 
the  coefficients  R  and  T  adopt  many  of  the  characteristics  of  the 
experimental  data.  In  any  event,  we  conduce  that  mere  refined  calcu¬ 
lations  should  use  a  more  realistic  excitation. 
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